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Abstract 

The dc conductance, gpj, through a finite Hubbard chain of size N (= 
1, 2, 3, . . .) connected to reservoirs is studied at T = in an electron-hole 
symmetric case. We calculate a spatial dependence of the self energy analyt- 
ically at u = within the second order in U, and obtain an inter-site Green's 
function Gjsrx, from which g^ can be determined, via the Dyson equation. 
The results depend strongly on whether ./V is even or odd. For odd N, a 
perfect transmission occurs, and gN = 2e 2 /h independent of the values of U. 
This may be attributed to a Kondo resonance appearing at uj = 0. For even 
N, gN decreases with increasing N, and converges to a finite constant which 
is a smooth decreasing function of U. These behaviors are essentially owing 
to the presence of the reservoirs, which make a quasi-particle description valid 
for low-energy states at uj <C tivp / '(Na); where vf is the Fermi velocity and 

a is the lattice constant. 

PACS numbers: 72.10.-d, 72.10.Bg, 73.40.-c 
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Effects of an electron-electron interaction U on the transport through a small system, 
such as a quantum dot and a quantum wire, have been a subject of current interest. So far, 
theoretical investigations for a quantum dot@'i and a quantum wirei seems to be developed 
separately because of the difference in their dimensionality. The purpose of this work is to 
study a system of intermediate size and investigate how the transport depends on the size. 
To this end, we consider a finite Hubbard chain of size N connected to Fermi-liquid reservoirs 
which are simulated by semi-infinite noninteracting leads (see Fig. |]), and calculate the dc 
conductance at T = with the perturbation approach! starting from a zero-dimensional 
limit. For small N, the system we consider may be regarded as a model for a quantum 
dot consisting of multi-levels, which has recently been examined by advanced numerical 
methods, such as a numerical renormalization group! and a quantum Monte Carlo method.! 
Specifically, in the present study, we concentrate on an electron-hole symmetric case. Thus, 
for large N, the system may be regarded as a model for a Mott-Hubbard insulator, which 
has been studied by a number of groups with the bosonization approachEHO] starting from a 
one-dimensional limit. 

We now proceed to details. The model consists of three regions (see Fig. [IJ); a finite 
Hubbard chain at sites 1 < % < N, and two semi-infinite leads on the left — oo < i < 0, and 
the right iV + l<i<+oo. The Hamiltonian is given by 

H = Hl + T~Cr + H° c + Hq 1 + H mix , (1) 
o 

K L = E [ ~ V S V ] C t*Cjcr , ( 2 ) 



i,j=—oo 
a 

oo 



= E [ ~ ^ 6 iJ ] C i* C jv ' ( 3 ) 

i,j=N+l 
a 

N r jj 



n° c = E 



-Uj + [ eo + 2" _ ) S v 



He* = U E 



N r 1 



1=1 



(5) 



Hmix = -V L E ( C L C 0a + H.C. ) - V R E ( Cn+1* C N<t + Hx - ) ■ ( 6 ) 



2 



Here c] CT (c jo .) is the creation (annihilation) operator for an electron with spin o at site j, 
rij a = CjpCjp, n is the chemical potential, e is the on-site energy. In each of the regions, 
the intra-region hopping matrix element — is finite only for the nearest-neighboring sites: 
—tij = —t [ Si + ij + ]. The connection between the interacting region and the left (right) 
lead is described by the tunneling matrix element Vl {vr). In the present study, we treat TCq* 
as a perturbation taking the unperturbed Hamiltonian to be TC^ = TCl + 'Hr + 1-L c + 7i. m i X . 
Especially, in a limiting case N — 1, our model corresponds to a single Anderson impurity, 
in which the perturbation expansion is valid for all values of [7.0 Our working hypothesis 
is that the ground state is changed continuously when N is increased to intermediate size 
owing to the presence of the noninteracting-leads.ill! 
We next consider the retarded Green's function; 

G js ,{u> + z0 + ) = -i fV ({c ja (t'), 4,(0)}) e^+^' , (7) 

<J 

where Cj a (t') = e lHt ' Cj a e~ lHt ' , (■ ■ ■) denotes the average Tr e _/3W ■ ■ • /Tre _/3W , (3 = 1/T, 
and the curly brackets denote the anticommutator. The spin index for the Green function 
has been omitted assuming the average to be independent of whether spin is up or down. 
Since the interaction is restricted for the electrons in the central region, the Dyson equation 
is written as 

N 

G^z) = G\f(z) + £ Gf(z) £^) G mj (z) , (8) 

l,m=l 

where Gf^ is the free Green's function corresponding to Ti,^°\ and £#(2) is the self-energy 
correction due to TVq 1 . Note that G{j = Gji and = S^j because of the time reversal 
symmetry of Ti. At T = 0, the dc conductance is expressed in terms of G^iioj + i0 + ) 



paw 



9n = ^- 4T L (0)r R (0) \G N1 (iO + )\ 2 , (9) 

where T Q (0) = kD(0)v* with a = L, R and D(0) = - Ji 2 / (2nt 2 ) being the local 
density of states at an edge of the isolated lead. The Green's function for u = 0, Gni(i0 + ), 
can be calculated solving the Dyson equation Eq. flD taking the self-energy E i:; (z0 + ) as 



an input. This part of the calculation can be regarded as a scattering problem of a free 
quasi-particle.0 It is described by an effective Hamiltonian H^J = Hl + 'Hr + T~ic + 
which has a renormalized matrix element, — tfj = —ty + Re£jj(0), in the central region 
Hq = Z)ij=i — tfj + ( e o + U/2 — y) 5ij c\ a Cj a . This mapping onto the free quasi-particle 
is justified under the condition lm'E i j(iO + ) = at T = 0, which can be shown with the 
perturbation theory .0 In other words, when the ground state is a Fermi-liquid, there exists a 
one-particle Hamiltonian which reproduces the Green's function exactly in the limit u = 0, 
T = 0. Thus, the physical quantities which are expressed in terms of GV,(iO + ), such as qn 
and the number of displaced electrons^ are reproduced exactly by using Ti$. 

We now calculate the self-energy. In what follows, we take the tunneling matrix el- 
ement to be vr = vl = t, for simplicity, and consider an electron-hole symmetric case, 
where the average number of electrons per site is unity, taking the parameters to be 
eo + U/2 — y = and /i = 0. In this case, the unperturbed Hamiltonian has a simple 
form = —t J2il-oo,cr ( c !+1ct c mt + H.c. V and effects of U starts from a second-order con- 
tribution illustrated in Fig. ^|. The corresponding retarded self-energy is written as 



V (J)/ , - n +N TT 2 [* dk 3 dk 2 dki 



i(ki+k 2 -k 3 )(j-j') 

(2nf 

x fkjiX ~ /fc 2 )(l ~ fkj) + (1 - fk 3 )fk 2 fk! QqN 

u + e k3 - e k2 - e kl + i0+ 
Here j and f are restricted in the central region 1 < < L, otherwise E^-)(a;-|-z0 + ) = 0, 
fk — + I] -1 , and e k = —2tcosk. It is straight forward to show that the imaginary 
part behaves as -ImEjy(w + i0 + ) oc uj 2 for small uj at T = 0. This behavior is quite 
different from that in a usual one-dimensional system, in which the imaginary part for a low- 
energy excitation at the Fermi point kp is proportional to |u;| as a result of the momentum 
conservation caused by the translational invariance. We next consider the real part at uj — 0, 
T = 0. Because of the electron-hole symmetry, the real part becomes zero when \j — j'\ 

(2) 

is an even number (= 2m); Re 2^(0) = 0. On the other hand, when \j — j \ is an odd 
number (= 2m + 1), Eq. ([Tt]) can be rewritten as; ReS 2m+ i (o) = -t {p(o)uys 2m+1 with 
p(0) = l/(2nt) being the density of states corresponding to HS Q \ and 



1 r' 2 j, t ; „ cos[A;i(2m + 1)] cos[A; 2 (2m+ 1)] cos[A;3(2m + 1)] 

S2m+i = — P / dk 3 dk 2 dk 1 — ■ — ■ 

2tt J-tt/2 cos/ci + cos/c 2 + cos/c 3 



Further, this expression can be simplified as 

(-l) m 7T 2 



S- 



2m+l 



cfeJ 2m+ l(z) 3{E 2m+ l(2)} -{J2m+l(^)}' 



(12) 



Here J2 m +i( z ) and E 2m+ i(;z) are the Bessel and Weber functionsrl The integration can be 
performed analytically by decomposing it into two parts)^ 



POO r- 

/ dz {J 2m +i{z)Y 
Jo 



1 3 1 3 1 

3_r 2 I -, —2m , 2m H — ; m, m H — ; — 

(2m+ 1)tt \2' 2' 2' 2 ' 2' 4 

where 3 F 2 (ai, a 2 , a 3 ; &i, fe 2 ; z) is the generalized hypergeometric functional and 



(13) 



r 2 2 

dz J2Tn+l(z) {J2m+l(z)} + {E 2m +l(z)} 



9 \ 2 m 

§) £<-i: 

/2Z\ 2 1 



// (2m + l) 2 
2/ \ / 



+ EE 

fc=l cj=l 



, (m + r 
[ 21 , 



+ 
21 



< 21 \ 2(2m+l) 
yl-k) [(2m-f 



l + (-D 



A; 2 



2(2m+ l) 2 


;(-i) fc +(-i) g ; 




[(2m+ l) 2 - (A; + g) 2 ; 


] [(2m + l) 2 - 


(A;-g) 2 )] 



(14) 



The result for S 2m +i is summarized in the Table 0, and is also plotted in Fig. [| The self- 
energy correction, 5 2m +i, is largest for the nearest-neighboring sites m = 0, and decreases 
rapidly with increasing m showing an oscillatory behavior. 

Using the above results for the self-energy, we determined the renormalized matrix ele- 
ment — tfj within the second order, and calculated GVi(iO + ) solving the scattering problem 
of a free quasi-particle described by HQ ■ Then, we obtained g^ from Eq. @. In Fig. f|, 

is plotted as a function of the size N for p(0)U = 1.0. When N is an odd number 
(= 2M + 1), the perfect transmission occurs, g 2 &i+i = 2e 2 /h, independent of the values 
of U. This remains valid even when the higher-order terms are included, which can be 
proved based on the quasi-particle description^ and the perfect transmission is understood 
as a result of the electron-hole and the inversion (vi = Vr) symmetries. Physically, it may 
be attributed to a Kondo resonance appearing at u = 0. In other words, this kind of an 



even-odd property is related to the level-structure of an isolated Hubbard chain of size N. 
For an odd N, the isolated chain has two eigenstates at u = 0, which corresponds to the 
Kramers doublet states. When the reservoirs are connected, the doublet states are changed 
to a Kondo resonance state and contribute to the perfect transmission. On the other hand, 
for an even N, there is no discrete eigenstate at uo = 0. These features reflect strongly on 
the transport of the connected system. 

When N is an even number, the dc conductance shown in Fig. § decreases with increasing 
N, and tends to a finite constant for large N. In Fig. || the dc conductance for the even 
case with N = 2M is plotted as a function of M for several values of p(0)U (= 0.5, 1.0, 1.5, 
and 2.0). The value of g 2 M converges well to a constant already at small M, and the 
constant decreases monotonically with increasing U. Figure ^| shows the U dependence of 
<?jv for iV = 50. Within the accuracy of the line thickness, the curve can be regarded as an 
extrapolated value owing to the early convergence. We note that the reduction of g 2 u from 
the universal value 2e 2 / h is the effects due to the electron correlation, and within the Hartree- 
Fock approximation g 2 M is unchanged from the universal value. In addition, our result shows 
a metallic behavior, i.e., g 2 M is finite in the limit M — > oo, in spite of the half-filled case. 
We have obtained this behavior by taking the limit T — > first, keeping M to be finite. In 
this limit, the sample and two leads are described by a single (ground state) wave function, 
and the phase coherence among the whole system plays an important role for the transport 
property. If the limit M — » oo is taken first at finite temperature, the phase coherence 
may be disturbed by thermal excitations. Qualitatively, our result is consistent with that 
obtained by Ponomarenko and Nagaosa,@ and by Odintsov, Tokura and Taruchai with the 
bosonization approach. Their results also show a metallic behavior bellow a characteristic 
temperature, and effects of a Mott-Hubbard gap on the dc conductance appear above the 
characteristic temperature. 

In the present study, we have studied a finite system connected to reservoirs described 
by the semi-infinite noninteracting leads. The presence of the reservoirs seems to make the 
quasi-particle description of a Fermi-liquid valid at lower temperatures than a characteristic 
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energy scale. For instance, in a limit N = 1, the system is reduced to a single Anderson 
impurity, and properties below the Kondo temperature are described by the Fermi-liquid 



tern of length L, a crossover seems to occur at a characteristic temperature To = Hvf/L 
(vp is the Fermi velocity) as pointed out by Kane and Fisher. i They have discussed that 
a Fermi-liquid behavior is expected at T <C To due to the cut off by the length L, and a 
Tomonaga-Luttinger behavior is expected at T 3> T>. In our case, the length is L = Na with 
a being the lattice constant. Since the characteristic temperature To is a decreasing func- 
tion of N, the low-temperature region described by a Fermi-liquid becomes small for large 
systems. Nevertheless, for small systems, the low-temperature region becomes relatively 
large. 

In summary, we have studied the dc conductance through a finite Hubbard chain con- 
nected to noninteracting leads based on the perturbation theory. The application of the 
quasi-particle description used in the present study is not restricted to the single-mode case, 
and will be applied to various systems which have a Fermi-liquid ground state. 

We would like to thank H. Ishii for valuable discussions. 




Analogously, when the reservoirs are connected to a finite one-dimensional sys- 
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TABLES 



TABLE I. Explicit expressions of S^m+i 



m S2m+1 

6 - \/3vr 

1 -14/5 + \/3vr/2 

2 290/147 - 5V3vr/14 

3 -164998/105105 + 2^/7 

4 2797322/2111655 - 22^3^/91 

5 -6632866/5731635 + llV3vr/52 
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FIGURES 



FIG. 1. Schematic picture of the model: (•) interacting region, (o) noninteracting leads. 



(2) 

FIG. 2. Feynman diagram for T}--,{z). 

FIG. 3. ReEg)(0) for odd (= 2m+l), normalized as ReS^ +1 (0) = -t {/9(0)f/} 2 S 2m +i 
with p(0) = l/(2vrt). The inset shows the region 2 < m < 11. Note that ReS^(O) = 0. 



FIG. 4. Size dependence of the conductance, <7at, for p(0)U = 1.0. 



FIG. 5. Conductance for even ./V (= 2M), g 2 M vs. M, for several values of p(0)U. 



FIG. 6. U dependence of for large even N (= 50). 
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